Viscous plasma evolution from gravity using AdS/CFT 
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We analyze the AdS/CFT dual geometry of an expanding boost-invariant plasma. We show 
that the requirement of nonsingularity of the dual geometry for leading and subasymptotic times 
predicts, without any further assumptions about gauge theory dynamics, hydrodynamic expansion 
of the plasma with viscosity coefficient exactly matching the one obtained earlier in the static case 
by Policastro, Son and Starinets. 
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One of the more challenging problems in theoretical physics is the understanding from first principles of the be- 
haviour of quark-gluon plasma, the phase of matter consisting of deconfmed quarks and gluons. There are strong 
indications that the plasma observed at RHIC is indeed strongly coupled (see e.g. [l[) and is well described by models 
based on hydrodynamics 0. This suggests the need for nonperturbative techniques to study the dynamics of the 
system. This need is especially acute if one would like to study non-equilibrium phenomena, thcrmalization etc. 

A very powerful technique for studying nonperturbative properties of gauge theory has emerged from string theory 
in recent years 0]. In its original form the AdS/CFT correspondence states the equivalence of Af = 4 Super Yang-Mills 
theory and string theory in a curved 10-dimensional AdS$ x S 5 background. What is crucial is that the string theory 
is easiest to handle in the regime of strong gauge theory coupling. 

Although the correspondence does not have a direct counterpart which works for QCD it has been argued that for 
studying features of the plasma not too far above the deconfinement phase transition it may be a good approximation 
as the plasma is strongly coupled and deconfmed. 

From a more general perspective, it is also interesting to study, for their own sake, dynamical time-dependent 
processes in the Af = 4 supersymmetric theory within the context of the AdS/CFT correspondence in order to have 
an example where such phenomenae can be calculated exactly, as well as to develop in this context new methods to 
address these dynamical issues. The methods used in the present paper translate questions on the behaviour of the 
plasma into certain questions within general relativity and we hope that this may inspire further research in both 
domains. 

Properties of the Af = 4 gauge theory at fixed finite temperature have been studied in some detail 0,@|, in particular 
shear viscosity has been calculated from perturbations around a static black hole background. On a more qualitative 
level, thermalization has been suggested to correspond to black hole formation in the dual description [6( , while cooling 
was advocated to correspond to black hole motion in the 5th direction [3] ■ 

In 0] a quantitative framework for studying such time-dependent phenomena has been proposed. The criterion of 
a nonsingular dual geometry was shown to pick out uniquely, in a boost-invariant setting, asymptotic perfect fluid 
hydrodynamical evolution for large proper-times. The resulting asymptotic geometry was shown to be analogous to 
a moving black hole. Further work within this framework includes [i^. fioj. 

The aim of this paper is to show that the criterion of nonsingularity predicts, when applied also to subasymptotic 
times, viscous hydrodynamic evolution with a specific viscosity coefficient. As a byproduct we obtain a nontrivial 
L* ■ consistency check of the AdS/CFT correspondence for the value of the shear viscosity. 
. £^ \ Boost-invariant viscous hydrodynamics. Let us consider the spacetime evolution of the energy-momentum 
tensor of an expanding fireball of plasma. The energy-momentum tensor is constrained by energy-momentum conser- 
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dfjT^" = (1) 
and, for the case of Af = 4 SYM theory that we consider here, tracelessness = 0. We will further restrict ourselves 



to boost-invariant evolution, first considered by Bjorken 11] as a model of the mid-rapidity region in heavy- ion 



collisions. This assumption is also commonly used in hydrodynamic simulations for RHIC j2J. We will further assume 
no dependence on transverse coordinates. It is natural to use the proper-time/spacetime rapidity coordinates for 
Minkowski space: 

ds 2 = -dr 2 + r 2 dy 2 + dx\ (2) 

Then energy-conservation and tracelessnes determine T^ v in terms of a single function - the energy density e(r) (for 
explicit expressions [IH see [§]). The dynamics of gauge theory should then determine e(r). In [§[ we proposed to 
use AdS/CFT to determine the proper-time dependence of energy density by first constructing the dual geometry to 
a given e(r) and then requiring its nonsingularity to fix the physical e(r). 



Let us review what would be the physical expectations for large proper-times if the gauge theory dynamics were 
described by viscous hydrodynamics. If there would be no viscosity, then we would have 



<t) = 4 (3) 
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If in addition we would have viscosity [l| 



V=^ (4) 

T 



then the energy density would behave like (see [Io| ) 



s r) = — T + . . . (5) 

with viscosity effects generating subleading deviations from the ideal fluid case ((3|). 

Construction of a dual AdS / CFT geometry. The procedure of constructing the dual geometry to a gauge 
theory configuration with given expectation value of the energy- momentum tensor was introduced in 12| . One adopts 
the Fefferman-Graham coordinates [l3T | for the 5-dimensional metric 

2 guudx^dx" + dz 2 
ds = - (6 

z A 

where the z coordinate is the 'fifth' coordinate while \x is a 4D index. One solves Einstein equations with negative 
cosmological constant: 

E a/ 3 = R a fj - -g a pR - 6g a/3 = (7) 
with a boundary condition for g M „ around z = 0: 

.. (8) 

The fourth order term is related to the expectation value of the energy-momentum tensor through 

N 2 

(Tfiv) = (9) 

The procedure is therefore to solve the 5-dimensional Einstein's equations (with negative cosmological constant A = 
—6) with the boundary condition (JSj> for a given spacetime profile of the gauge-theoretical (!),„). In the following since 
we will be dealing directly with the metric and hence with gffi so we will supress the factor N 2 / (2tt 2 ) throughout the 
computation reinstating it only in the final discussion. 

Nonsingularity and viscosity. The metric consistent with the symmetries of the boost invariant expansion has 
the form 



-e 



a{r,z) dT 2 + T 2 e b(r,z) dy 2 + e c(r.z) dx 2^ + dz 
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ds z = — : — ^ (10) 

Z A 

In Q the dual geometry was determined for asymptotic times with the energy-density behaving like 

e(r) = -4 (11) 

for large proper-times. The resulting metric coefficients were given for large proper-times as functions of the scaling 
variable 

«=4 (12) 

T 4 

It was found in [1] that the resulting geometry was nonsingular only when s = 4/3, thus corresponding to perfect 
fluid hydrodynamics. Of course subleading corrections like §5§ are possible. The explicit leading coefficients for this 
case were found to be 

(l-^ 4 /3) 2 
a{v) = lQg 1 + .V3 

b(v) = c(v) = log(l+w 4 /3) (13) 



The resulting geometry looks like a black hole whose horizon (in the Fefferman- Graham coordinates) moves in the 
fifth dimension as zq = 3~T3. This leads to the temperature 

r=^I = 4r-* (14) 

7T Zq 7T34 

The coefficient \/2 comes from the special form of the black hole metric in Fefferman-Graham coordinates. See § for 
a discussion. 

In an interesting recent paper Nakamura and Sin [Io| determined the leading, subasymptotic in proper-time, 
corrections to the metric coefficients like 

a(r, z) ~ a(v) + ai(v)— + . . . (15) 

T 3 

which represent e(r) of the form (|5|). The coefficients which can be extracted from the results of in a form 
convenient for proceeding to higher orders are 

(9 + t ,> 4 
ai{v) = 2m 9~v» 

v 4 3 — u 4 

ci(v) = -2?7o— — j - 77 log— (16) 
3 + v 4 3 + w 4 

However to this order in r (0(t~9)) the Riemann tensor squared 

m 2 = R^ al3 R^ aP (17) 
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was found to be finite for any t]q, the first divergence appearing only at order r~a . This behaviour suggested that in 
order to determine the coefficient of viscosity 770 one had to go one order higher in t~3 . Namely one has to find the 
metric coefficients to second order like 

a(r, z) = a(v) + ax{v)\- + a 2 (v)-r + • • • (18) 
r 3 r 3 

In order to setup a systematic expansion procedure one expands the left hand side of the Einstein equations Ea = 

2 4 2 4 2 '- 1 

(•7-3 E tt ,t3 E tz ,t~3 E zz ,t~s E yy ,T3 E xx ) in powers of t~S: 

E A = Ef(v) + 4%) 4 + 4 2) («)4 + • • • (19) 

T 3 T3 

where A = 1 ... 5 numbers the five nontrivial components of the Einstein equations mentioned above. The prefactors 
are chosen in such a way as to have a uniform expansion. The leading order results (| 13|) satisfy Ejj^ — while 

the first subleading corrections (fTB]) satisfy E^ = 0. Therefore in order to find 1x2 (v), ^(v), C2(v) we have to solve 
(2) 

the equations iv^ = with the boundary condition that these coefficient functions should vanish at v = 0. The 
resulting expressions depend on the viscosity coefficient 770 and another coefficient which we denote by C. The explicit 
expressions are 

r /4w 5 (27 + 9u; 4 + 2u; 8 ) „ „, w 3 (405 + 171w 4 + 189w 8 + 5w 12 - 2w w ) „ w 3 (9 + 2w 4 + w s ) \ , 
^ v) = J { 3(9 - W 8 ) 2 ~ 24 ^— (9^ ° O-^ 8 ) 2 
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&a(«) = -2ca(t>) + — -r — ^ + r? 2 ; Q ■ ,, 2 + -t? 2 log — — + C- 
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6 + 2w 4 2\/3 ° (3 + « 4 ) 2 2 " J b 3 + « 4 12(3 + v 4 ) 

T / 4w 9 (9 + w 4 ) 8w 3 arctanh^ 2 24w 3 (w 4 - 15)(3 + 5w 4 ) w 3 (l + w 4 ) , , 
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One can express these integrals in terms of elementary functions and dilogarithms but the expressions are rather 
lengthy and will not be presented here. 



It remains to determine when the background geometry given by the coefficients up to second order is nonsingular. 
To this end we calculate the curvature invariant defined by (|17p and expand it in the scaling limit up to the order 

4 

t~3 . The resulting expression at this order has the form 

2 . 1 polynomial in v, t]q and C , . 

m = nonsingular terms + — — ^ (21) 

T3 (o — V 4 ) 4 (3 + V 4 ) 

We see that there is a potential singularity at v — 3^. It turns out that the singularity is cancelled exactly when 

% 2 - ^ (22) 
There is no restriction on C at this order. We expect one would have to perform the analysis to the next order to fix 



C [14|. 



The above calculation shows that nonsingular dual geometry is not possible for an exact perfect fluid, but that 
viscosity effects in the proper-time evolution are present with a uniquely fixed value of the viscosity coefficient given 
by (f22|) . It is interesting to compare this value with the shear viscosity obtained by Policastro, Son and Starinets Q 
who derived it by studying the reponse of a static plasma at fixed temperature to small perturbations. 

To this end let us take the value of (shear) viscosity obtained in Q at a fixed temperature T: 

r, = — a = -N 2 C T 3 (23) 

If we insert the proper-time dependence of the temperature (|14l) for the evolving plasma into the above expression we 
obtain 

N c 11 

where we have factored out the coefficient N 2 /(2ir 2 ) appearing in ((9J. The resulting estimate for the viscosity 
coefficient r/Q is therefore 

^o = ^3-|^| (25) 

which is exactly the value (|22p for which the dual background geometry of the evolving plasma is nonsingular. 

Discussion. In this paper we have studied the proper-time evolution of a boost-invariant plasma using the 
AdS/CFT correspondence. We have shown that the requirement that the dual geometry is nonsingular predicts the 
proper-time evolution of the energy density to be equal [ItJ to the one found from viscous hydrodynamics with the 
viscosity being exactly the one following from rj/s = l/(47r) in the static case. The computation involves the nonlinear 
regime of gravity on the AdS/CFT side and it is encouraging for other possible applications that the method can 
capture such fine details of the gauge theory dynamics. 

It would be very interesting to study in more detail the features of this geometry and its thermodynamics as 
well as to apply these techniques to other dynamical non-equilibrium processes. 
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